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A ballistic strip of graphene (width W 3> length L) connecting two normal metal contacts is known 
to have a minimum conductivity of 4e 2 /irh at the Dirac point of charge neutrality. We calculate 
what happens if one of the two contacts becomes superconducting. While the ballistic conductance 
away from the Dirac point is increased by Andreev reflection at the normal-superconductor (NS) 
interface, we find that the minimum conductivity stays the same. This is explained as a manifestation 
of pseudo-diffusive conduction at the Dirac point. As a generalization of our results for a ballistic 
system, we provide a relation between the conductance Gns of an arbitrarily disordered normal- 
superconductor junction in graphene and its value Gn when both contacts are in the normal state. 

PACS numbers: 74.45.+C, 73.23.-b, 74.50.+r, 74.78.Na 



I. INTRODUCTION 

The effect of a superconducting contact on the conduc- 
tance of a metal is qualitatively different if the electrons 
propagate ballistically or diffusively through the metal: 
A superconducting contact doubles the conductance of a 
ballistic metal, while the conductance of a diffusive metal 
remains the samei The conductance doubling in a bal- 
listic metal happens because an electron incident on the 
superconductor is Andreev reflected as a hole of oppo- 
site charge, hence doubling the current. The effective 
length L of the conductor is doubled as well, because the 
hole has to make its way back, and this effect cancels the 
conductance doubling in the case of a diffusive metal. 

These two competing effects are summarized by the 
approximate relation 

Gns (L) « 2G n (2L), (1) 

between the conductance Gns of a normal-metal- 
superconductor junction and its normal-state value Gn- 
As derived in Ref. 2 from a generalized Landauer for- 
mula, the relation holds if quantum corrections of 
order e 2 /h can be neglected. It also requires an ideal 
NS interface, without a mismatch of Fermi wave lengths 
between N and S. 

A novel transport regime called "pseudo-diffusive" ap- 
pears in an impurity-free carbon monolayer (= graphene) 
at zero carrier concentration (= Dirac point). Because of 
the vanishing density of states the transmission through 
a strip of undoped graphene (width W, length L) oc- 
curs entirely via evanescent modes. For a short and 
wide strip there is a large number N e g = W/L ^> 1 
of evanescent modes with transmission probability of or- 
der unity (open channels). Several recent investigations 
have found that the transport properties of undoped bal- 
listic graphene are similar to those of a diffusive metal 
with the same number N e g of open channels: Both sys- 
tems have the same conductance 3,4 G ~ N c ge 2 /h and 
the same shot noise powes* (Fano factor 1/3). In a 
Josephson junction both systems have the same current- 
phase relationship, critical current, and current-voltage 



characteristic!^ This correspondence between evanes- 
cent modes in graphene and diffusion modes in a dis- 
ordered metal is not limited to a carbon monolayer, but 
applies to a bilayer as wel\X The effect of disorder on the 
evanescent modes has also been studied^ 

In this paper we present one more manifestation of 
pseudo-diffusion in graphene, by calculating the ratio 
Gns/Gn- Far from the Dirac point, at high carrier con- 
centrations, this ratio is between 1 and 2 — as expected 
for a ballistic metal with a Fermi wave length mismatch. 
At the Dirac point, however, the ratio approaches unity 
— as in a diffusive metal, but without any scattering by 
impurities or lattice defects. 

The outline of this paper is as follows. In Sec.^lwe de- 
scribe the scattering problem of a lightly doped graphene 
strip between two heavily doped metal contacts. In Sec. 
IIIII we calculate the conductance Gns of the NS junction 
from the probability of Andreev reflection, as a function 
of the Fermi energy in the graphene strip. We compare 
with the known results 3 ^ for Gn in Sec. lIVI demonstrat- 
ing the equality Gns = Gn at the Dirac point. This 
explicit calculation is for a ballistic conductor. In the 
Appendix we derive a more general relation between Gns 
and Gn, valid also for a disordered conductor. 



II. SCATTERING MATRIX 

We consider a graphene sheet having width W in the 
y-direction, with regions x < (S region) and x > L 
(N region) covered by superconducting and normal elec- 
trodes, respectively. The central region < x < L is 
lightly doped, with Fermi energy /i relative to the Dirac 
point. The contact regions x < and x > L are heav- 
ily doped, with Fermi energy fi^ much larger than both 
fi and the superconducting gap Aq. An electron inci- 
dent on the superconductor at an energy e just above 
the Fermi level is Andreev reflected as a hole at an en- 
ergy e below the Fermi level. (See Fig. [I]) The re- 
flection occurs predominantly via propagating modes for 
\x 3> hv/ min(VF, L) (with v the energy-independent car- 



2 




FIG. 1: Top panel: geometry of a graphene sheet (G) with 
one normal (N) and one superconducting contact (S). Bottom 
panel: Spatial dependence of the Dirac point (solid line). The 
Dirac point is far below the Fermi level (at zero energy) in 
the heavily doped contact regions x < and x > L. In the 
lightly doped central region < x < L the Dirac point is 
much closer to the Fermi level. An electron incident from 
the normal contact (N) above the Fermi level (filled circle) 
is Andreev reflected by the superconducting contact (S) as a 
hole below the Fermi level (open circle). In the limit /«, s — ► 
the electrons and holes are transmitted through the central 
region via evanescent modes. 



rier velocity in graphene). That regime was studied in 
Refs. IgllCl In the limit \i — ► at fixed W, L only evanes- 
cent modes remain. This is the regime studied here. 

We calculate the scattering matrix of the NS junction 
using the Dirac-Bogoliubov-de-Gennes equation^ 



vp ■ cr + U 

A* 



A 

-vp ■ cr 



U 



(2) 



where p — —ih(d/dx,d/dy) is the momentum operator 
and cr = (a x ,cry) is the isospin operator. The excitation 
energy e > is measured relative to the Fermi level (set 
at zero). The electrostatic potential U and pair potential 
A have step function profiles, 



U(x) = 
A(x) = 



—fioo if x < or x > L, 
—fi if < x < L, 

A if x < 0, 
if x > 0. 



(3) 
(4) 



The Fermi wave length = hv/fioo in the contact 
regions is sent to zero by taking the limit fioo — > oo. 
The infinitely large mismatch with the Fermi wave length 
Xf = hv/n in the central region would fully suppress An- 
dreev reflection of nonrelativistic electronsiii It is a spe- 
cial property of the relativistic wave equation J3J) that a 
nonzero Andreev reflection probability remains regard- 
less of the Fermi wave length mismatch^ 

For aspect ratios W ^> L the boundary conditions in 
the y-direction are irrelevant. We take periodic bound- 
ary conditions for simplicity. Different wave vectors 
q n = 2im/W (with n = 0, ±1, ±2, . . .) in the y-direction 
are not coupled, so we can consider each transverse mode 
separately. In the normal regions there are four eigen- 
states of Eq. (J2J for each n and e, corresponding to elec- 
trons and holes propagating in either +x or — x direction. 
In the superconducting region there are two eigenstates 
that decay for x — * — oo. These eigenstates are given in 
Ref. 0. By matching them at x — and x — L for an 
electron incident from the normal contact we obtain the 
reflection amplitudes Tee{s) and rj^(e). Matching for 

an incident hole gives the amplitudes rjj^ (e) and r^(e). 
Together these four reflection amplitudes determine the 
scattering matrix 



The result of this calculation is 
1 



Teh 



The 



X 



cos a e cos ah, 



(5) 



(6a) 
(6b) 



r ee = YJX, r hh = Y h /X, 
X = e 2l/3 (cos a e cos # e — z sin6 l e )(cos a^, cos 6*^ — ismOh) 

— sin a e sin 9 e sin ah sin 6h , (6c) 
Y e = e 2 ^ sin a e sin 9 e (cos ah cos 6h — i sin Oh ) 

— sine*/! sin Oh (cos a e cos e + i smO e ), (6d) 
Yh = e 2l>3 smahSmOh{cosa e cosO e — ismO e ) 

— sin a e sin e (cos ah cos Oh + i sin Oh ) ■ (6e) 



We have defined the angles 
hvq. 



/' 



at, 



. ( kvq n 
arcsm f 



\e — ji 



.(7a) 



L{£ + /-*) „ Lie - u) .„ . 

cosa e , Oh = — cosa/j, (7b) 



fiv 

(3 = arccos(£/Ao) 



tiv 



(7c) 



One can verify that S^ n '(e) is a unitary matrix for e < 
A . 



III. CONDUCTANCE 



Because there is only a single superconductor, we can 
take Aq real without loss of generality. 



The linear response conductance Gns = 
limy^o dl/dV is given by the Blonder-Tinkham- 
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FIG. 2: Conductance of the NS junction versus Fermi energy, 
calculated from Eqs. JUJ and I1U1 . The dotted lines are the 
asymptotes (I1H . 



' he 



(0)| 



Klapwijk formula 12 

Gns = go^-lrZHW + lrfcW 

n 

= 2 5 o£h ( " ) < 



(8) 



(In the second equality we used the unitarity of the scat- 
tering matrix.) The conductance quantum is go = 4e 2 //i, 
where the factor of 4 accounts for the twofold spin and 
valley degeneracies in graphene. The sum over modes 
runs from — fx^W/hv to +/i QO W/hv. In the regime 
fioo — > oo, W/L — » oo of interest here, the sum over 
modes may be replaced by an integration over the trans- 
verse wave vector q, 



f°° W 
Gns = 2g / Ra(q) tt - dq. 

J-oc 27T 



(9) 



with Rji{q n ) = \r he (0)\ the probability for Andreev re- 
flection at the Fermi level. 

From Eq. @ we obtain the expression 



[{n/hv) 2 - q 2 cos(2fcL)]' 



(10) 



yj (n/hv) 2 — q 2 the longitudinal wave vector. 
> \fJ.\/hv the longitudinal wave vector is imag- 

Substitu- 



with k 
For \q\ 

inary, corresponding to an evanescent mode, 
tion of Eq. ljTU|l into Eq. gives the conductance versus 
Fermi energy plotted in Fig. [2J The asymptotes are 



G 



NS 



TT^goW/L 
0.38g \n\W/hi 



for 
for 



< hv/L. 
> hv/L. 



(11) 



IV. COMPARISON WITH A NORMAL 
JUNCTION 

In Refs. and the geometry of Fig. ^ wa s studied 
in the case that both contact regions are in the normal 
state. The transmission probability T from one contact 
to the other in that case was found to be given bj«^ 



T(q) 



k 2 cos 2 (kL) + {n/hv) 2 sm 2 (kL) 



(12) 



As before, q is the transverse wave vector and k = 
\J (ji/hv) 2 — q 2 the longitudinal wave vector. In the limit 
fj, — > of undoped graphene we have T = 1/ cosh 2 (gi), 
while Ra = 1/ cosh 2 (2gi). The effective doubling of the 
length in the case of a superconducting contact is consis- 
tent with the physical picture described in the introduc- 
tion. 

For arbitrary /z, we find by comparing Eqs. i|l(J[l and 
(|12fl that the probability for Andreev reflection Ra is 
related to the normal state transmission probability T 
by 



Ra = 



T 2 



(2-T) 2 



(13) 



This relation, derived here from the Dirac-Bogoliubov- 
De Gennes equation, is the same relation as for the usual 
(nonrelativistic) Bogoliubov-De Gennes equation. 2 In the 
Appendix we show that this relation is not restricted to 
ballistic NS junctions, but can be generalized to arbitrary 
disorder in the normal region. 
The conductance 



r°° w 

Gn = 50 / T(q) — dq 



in the normal case has asymptote: 



3,4 



G 



'N 



-'goW/L 
:go\n\W/hv 



for 
for 



< hv/L, 
> hv/L. 



(14) 



(15) 



Comparison with the asymptotes l|llfl shows that the 
conductance is enhanced at large when one of the 
contacts becomes superconducting, as expected for a bal- 
listic junction. The enhancement by a factor 1.5 is be- 
low the factor of two enhancement expected for an ideal 
NS interface, because of the mismatch of Fermi wave 
lengths. For /i — > the enhancement disappears and 
Gns becomes precisely equal to Gn, as a manifestation 
of pseudo-diffusive conduction. 

Since the first experiments on induced superconduc- 
tivity and Andreev reflection in graphene have now been 
reported^*^ our theoretical predictions may well be 
tested in the near future. 
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APPENDIX A: GENERAL RELATION 
BETWEEN THE CONDUCTANCE OF AN NS 
JUNCTION IN GRAPHENE AND ITS 
NORMAL-STATE VALUE 

In Ref. a general relation was derived between the 
conductance Gns of an NS junction and its normal-state 
value Gn, on the basis of the Bogoliubov-De Gcnncs 
(BdG) equation. 15 Here we show that the same relation 
holds also for an NS junction in graphene, described by 
the Dirac-Bogoliubov-De Gennes (DBdG) equation^ 

The relation of Ref. is expressed by the pair of equa- 
tions 

G NS = 2g ^ ^ G N =.go]T T «' ( A1 ) 

n ^ n ' n 

where go — se 2 /h is the conductance quantum (including 
an s-fold spin and valley degeneracy factor). The trans- 
mission eigenvalues T n are the eigenvalues of the matrix 
product ttf , with t the transmission matrix at the Fermi 
level when both contacts are in the normal state. The 
region in between the two contacts may be arbitrarily 
disordered, so t is no longer a diagonal matrix [as in the 
ballistic case considered in the text, leading to Eq. lfT3jl], 
We also allow for intervalley scattering, in which case go 
includes only the spin degeneracy (s = 2). 

The relation (|A1|) relies on a separation of length scales 
on the superconducting side of the NS interface: The 
first length scale = hv/noa determines the distance 
over which the step Q in the electrostatic potential U (x) 
affects the wave functions in S. The second length scale 
£ = hv/Ao determines the distance over which the step 
(0J in the pair potential A(x) affects the wave functions 
in S. Since A^ <C £, we may treat the scattering by U 
and A independently of each other. 

Setting first A = 0, we consider the normal-state scat- 
tering problem, including the electrostatic potential U 
and any disorder in the region x > 0. This determines 
the scattering matrix Sq{s) for electrons and its counter- 



part Sq(— e) for holes. There is no mixing of electrons 
and holes for A = 0, so the full normal-state scattering 
matrix has the form 

Ssie) = [ S f s ^_ £) ) • (A2) 

The basis of hole states is chosen such that it is 
the time reversed of the basis of electron states 
^h{e) = (c z C§) r x )4'*(— s), with Pauli matrices <Ji and n 
acting, respectively, on the isospin and valley degree of 
freedom^ The normal-state reflection and transmission 
matrices are submatrices of So, 

So = . (A3) 

Since we assume zero magnetic field, So is symmetric as 
well as unitary. 

Setting then U = 0, we consider the Andreev reflection 
by the superconductor at x < (in the absence of any 
disorder for x > 0). The scattering matrix for Andreev 
reflection at subgap energies e < Ao is given by 

S A (s) = C(e) ' C(e) = exp[-iarccos(e/A )]. 

(A4) 

This has the same form in the BdG and DBdG equation. 

The total reflection matrix S is constructed from Sn 
and Sa in the same way as in Ref. y- We need the sub- 
matrix 

r he (e) = C(e)t*(s)[l - C 2 (e)r(e)r*(- £ )]- 1 <'( £ ). (A5) 

Setting e = and using the symmetry and unitarity of 
So, the conductances 

G NS - 2 9o Trr he (0)rj; e (0), G N = g Tr t(0)t+ (0) (A6) 
take the form of Eq. IjAip . 



1 Two reviews of conduction through normal-metal- 
superconductor junctions are: H. Courtois, P. Charlat, Ph. 
Gandit, D. Mailly, and B. Pannetier, J. Low Temp. Phys. 
116, 187 (1999); C. W. J. Beenakker, Rev. Mod. Phys. 
69, 731 (1997). The first review is from the perspective 
of semiclassical theory, the second from the perspective of 
random-matrix theory. 

2 C. W. J. Beenakker, Phys. Rev. B 46, 12841 (1992). 

3 M. I. Katsnelson, Eur. Phys. J. B 51, 157 (2006). 

4 J. Tworzydlo, B. Trauzettel, M. Titov, A. Rycerz, and C. 
W. J. Beenakker, Phys. Rev. Lett. 96, 246802 (2006). 

5 M. Titov and C. W. J. Beenakker, Phys. Rev. B 74, 
041401 (R) (2006). 

6 J. C. Cuevas and A. Levy Yeyati, |cond~m at/0609138 

7 I. Snyman and C. W. J. Beenakker, |cond -mat /0609243 



J. A. Verges , F. Guinea, G. Chiappe, and E. Louis, 
cond-mat/0610201 
9 C. W. J. Beenakker, Phys. Rev. Lett. 97, 067007 
(2006). The eigenstates are given in the appendix of 
cond-mat/0604594 

10 S. Bhattacharjee and K. Sengupta, cond-mat/0607489 

11 A. F. Andreev, Sov. Phys. JETP 19, 1228 (1964). 

12 G. E. Blonder, M. Tinkham, and T. M. Klapwijk, Phys. 
Rev. B 25, 4515 (1982). 

13 H. Heersche, P. Jarillo-Herrero, J. B. Oostinga, L. M. K. 
Vandersypen, and A. F. Morpurgo, preprint. 

14 A. Geim, private communication. 

15 P. G. de Gennes, Superconductivity of Metals and Alloys 
(Benjamin, New York, 1966). 



